Evolution & Development

Regenerative capacity varies widely across flatworms (Platyhelminthes). Whereas catenulids, microstomids and planarians can regenerate a complete head de novo, other flatworms cannot. This striking diversity raises a longstanding evolutionary question: does whole-body regeneration represent an ancestral trait that was subsequently lost in multiple lineages, or did it evolve convergently? Addressing this question requires comparative analyses of the molecular mechanisms underlying regeneration across phylogenetically diverse flatworms. Here, we focus on Wnt signaling, a deeply conserved regulator of antero-posterior (A-P) patterning and a central determinant of head-versus-tail identity during planarian regeneration, to establish a mechanistic framework for such comparisons. Although Wnt signaling has been studied extensively in planarians and parasitic neodermatans, its evolution and deployment in other flatworm clades remain poorly characterized. To address this gap, we characterized the complement of Wnt signaling components in two early-diverging flatworm clades, Catenulida and Macrostomorpha, with particular emphasis on expression and function in the catenulid Stenostomum brevipharyngium. Phylogenetic analyses reveal the ancient loss of six Wnt families and one secreted Frizzled-related protein (sFRP) family in the last common ancestor of flatworms, followed by additional lineage-specific gene losses and expansions. Moreover, several Wnt pathway components display markedly divergent expression patterns between catenulids and other flatworms, while functional analyses indicate corresponding differences in their regenerative deployment. Together, our findings reveal a dynamic evolutionary history of the flatworm Wnt signaling toolkit and establish a comparative framework for testing whether the molecular circuitry underlying head regeneration is ancestrally conserved or has evolved independently in distinct flatworm lineages.

Sea urchins are pivotal models in cellular and developmental biology, but their biphasic life cycle with an extended larval life limits the ability to study post-metamorphosis and adult characters. Here, we introduce a genomically-enabled model system, the tuxedo urchin Mespilia globulus, which has rapid access to late life stages in sea urchins. We describe how we cultured M. globulus in a landlocked marine facility, raised larvae under artificial conditions and closed their life cycle. We established the experimental tractability of M. globulus: we labelled transcripts by hybridization chain reaction (HCR), and knocked out pigmentation genes to produce albino larvae using CRISPR/Cas9. We generated chromosome-scale genome assemblies for two individuals representing both sexes and two color morphs (red and blue), and compared the organisation of the 21 chromosomes of M. globulus with that of other camarodont echinoid models. We determined that M. globulus showed a conservative gene repertoire lacking the gene family expansions seen in other camarodont sea urchins. We annotated the complement of genes associated with pigmentation, immune and nervous systems and profiled their expressions in tissues and organs. Finally, we surveyed sex-related regions in genomes using genome assemblies and resequencing data, finding no evidence of heteromorphic sex chromosomes in M. globulus. Our findings highlight the accessibility of this new sea urchin model for studying the metamorphosis and adult biology of sea urchins.

Sexual dimorphism is a well-studied phenomenon in which males and females of a species exhibit different phenotypic characteristics. However, the degree of sexual dimorphism within a species may vary among populations and may be influenced by geographic factors such as seasonality, prey availability, and population density. In this study, we examined the effects of sexual dimorphism and geographic region on the craniomandibular morphology of bobcats (Lynx rufus) in Washington State. Using digital calipers, we measured 13 craniomandibular traits from 35 males and 42 females. We found that the degree of skull dimorphism does not differ between the eastern and western side of the Cascade Mountains. However, we found significant effects of sex and geographic region on craniomandibular morphology. Across the geographic regions, males were significantly larger than females in most craniomandibular traits, excluding postorbital constriction breadth (POC). After correcting for expected size differences, only condylobasal length (CBL) was disproportionately longer in males, indicating subtle shape dimorphism within this trait, that is largely unexplained by allometry. Mechanical advantage (MA) of the temporalis and masseter muscles did not differ between the sexes at the canine or the molar bite points. Comparisons between populations east and west of the Cascade Range revealed that while sex remains the primary indicator of overall size, there is a small regional effect on size as western bobcats are slightly larger than eastern bobcats. These regional differences highlight how environmental pressures including local prey availability and severe winter conditions may influence functional skull design. Future ecological and behavioral studies are needed to further clarify how these localized environmental drivers shape phenotypic variation across populations.

Hard cornification is the developmental process that produces the hardened epidermal structures that characterize tetrapod diversity, including scales, feathers, hair, and claws. This process relies on keratin proteins and additional reinforcing proteins that assemble, cross-link, and strengthen tissues to become mechanically resilient. Although the molecular composition of hard cornified tissues is known to vary among tetrapods, the program is least understood in amphibians. Within frogs, hard cornification occurs in the mouthparts of most tadpoles and in the nuptial pads, hindfoot claws, and foot tubercles of some adults. We characterized transcriptomic profiles of Bombina orientalis tadpoles at three developmental stages before, during, and after mouthpart cornification. We identified 25 keratin genes orthologous to amniote hard keratins in the B. bombina genome, 18 of which are upregulated in cornifying larval mouthparts. Furthermore, two of these hard keratin paralogs and an amphibian epidermal differentiation gene exhibit distinct spatial expression patterns within the jaw sheaths and keratodonts evidenced by in situ hybridization, suggesting they have discrete functional roles. The newly discovered diversity of the hard keratin repertoire in Bombina , and the concurrent expression of many of these paralogs in tadpole mouthparts, suggests that anuran hard cornification is much more complex than adult cornified traits thus far characterized in Xenopus frogs. Our results also challenge the notion that the diversification of keratins in Amphibia is primarily associated with the transition to more terrestrial habits, and may instead be shaped, in part, by the novel feeding apparatus of aquatic tadpoles.

Color polymorphisms offer a powerful lens for examining how selection shapes genetic and phenotypic variation. Genomic analyses of color variation provide deeper insight into how genetic differences and phenotypic plasticity contribute to trait variation. Here, we use the bluefin killifish (Lucania goodei) to study this relationship. Male fins can display either red or yellow coloration based on pigments (pterins) or blue structural coloration (iridophores) often induced in UV-depauperate environments. Female fins lack pigmentation. To characterize the molecular basis of this variation, we assembled a high-quality genome and examined gene expression in the anal fins of male and female bluefin killifish to identify the molecular mechanisms underlying color variation. Our gene expression analysis found sex differences driven by male-biased upregulation of melanin synthesis and canonical pigment pathway genes. We additionally identified that the progesterone sex receptor was upregulated in females. Among male color morphs, several genes involved in xanthophore and iridophore formation differed between blue and non-blue males, including the transcription factor tfec, and iridophore-patterning genes atic, pnp4a, and sox10. Extraocular opsins, which are involved in non-image-forming light perception, were expressed at low levels across all fins, but showed population- and morph-specific patterns; notably, rgrb (retinal G protein coupled receptor b) was upregulated in blue males. To the best of our knowledge, this is the first report of an extraocular opsin varying as a function of population and color morph. Together, these results identify the pigment pathways underlying male-limited color polymorphisms and reveal candidate genes likely involved in light-sensitive phenotypic plasticity.

The red palm weevil, Rhynchophorus ferrugineus, is the most economically destructive palm pest worldwide, threatening livelihoods, food security, and ecosystems across 49 countries. Weevil management currently relies predominantly on chemical insecticides, raising significant environmental and public health concerns. Despite its global agricultural importance, genetic approaches to pest management and the mechanistic basis of genome-editing strategies in Rhynchophorus remain largely unexplored. Here, we employed CRISPR/Cas9 genome editing to disrupt the R. ferrugineus ommochrome biosynthetic pathway -- a multi-enzymatic metabolic cascade that converts tryptophan into ommochrome pigments, including brown, yellow, and red pigments. We targeted two key pathway components: the ATP-binding cassette transporter white and the heme peroxidase cardinal. Both genes were ubiquitously expressed, with peak expression levels in the gut, fat body, and head. Elevated transcript levels were observed across early, mid, and late pupal stages and in 0-, 1-, and 2-day-old adult males and females, consistent with the progression of eye pigmentation throughout the R. ferrugineus life cycle. Embryonic microinjection of a single guide RNA (sgRNA)-Cas9 ribonucleoprotein complex targeting white produced in the Generation-0 (G0) adults with a distinct, white-eyed phenotype with a brownish outer margin, in contrast to the black eyes of wild-type adults. Genome-edited cardinal mutant adults displayed a translucent, brownish-white-eyed phenotype, with white streaks that gradually transitioned to a persistent translucent reddish-brown eye coloration. Mutations in both genes were confirmed in G0 adults by genomic DNA sequencing. Mutant adults were crossed to generate heterozygous G1 (+/-), G2 (-/-, -/+, and +/+), and G3 lines (-/-) with genotypes verified as carrying 2-, 3-, 9-, and 13-nucleotide deletions. A stable, heritable eye-color phenotype was established in homozygous knockout (-/-) G3 lines for both white and cardinal, confirmed by unambiguous indel (insertions/deletions) genotyping. Inheritance analysis revealed that both genes are X-linked, following a classical Mendelian sex-linked pattern: paternal alleles are transmitted exclusively to daughters, while maternal alleles are inherited equally by both daughters and sons. This study establishes the first fully homozygous knockout strain in R. ferrugineus and, by characterizing sex-linked inheritance in a coleopteran system, advances our understanding of how CRISPR/Cas9 can be efficiently applied to destructive palm weevil species. The present study represents the first report of CRISPR/Cas9 genome editing in any weevil (Curculionidae), using white and cardinal as marker genes. These findings provide a valuable platform for functional genomics and genome engineering in R. ferrugineus and offer a translational framework for genome editing in the invasive South American palm weevil, R. palmarum, laying a solid foundation for the development of gene-drive strategies aimed at sustainable palm weevil population control.

Quantifying skeletal mineralization phenotypes in larval fish is complicated by the natural curvature of the notochord and by sample-to-sample variability in orientation, staining and imaging. Consequently, many studies rely on summary measures such as vertebral counts or total stain intensity. Here we present SCAMP (Spinal Calcification & Mineralization Profiler), an open-source, GUI-based Python tool that computationally straightens the curved notochord of Alizarin Red S-stained fish larvae and generates standardized mineralization profiles along the spinal axis. This approach reduces positional and shape variability, allowing direct, quantitative comparison of calcification patterns within and between experimental cohorts, without requiring programming expertise. We validate SCAMP using a zebrafish model of Pseudoxanthoma elasticum (abcc6aelu15/elu15), recovering genotype-specific differences in the intensity, extent and spatial distribution of ectopic calcification. Using SCAMP, we further show that inorganic pyrophosphate (PPi) supplementation of the medium suppresses ectopic notochord calcification, alters the anterior-posterior distribution of mineralized regions in homozygous mutants, and promotes mineralization at physiological vertebral sites. We also show that methylene blue, a routine antifungal additive in fish medium, reduces baseline calcification, with the most pronounced effects observed in heterozygous controls. SCAMP is freely available and has the potential to be adapted to other fish species used in skeletal and mineralization research.

Hox genes play central roles in specifying positional identities along the vertebrate anterior-posterior axis. In mice, genetic analyses have demonstrated that Hox genes distributed among the four Hox clusters contribute to vertebral patterning, with extensive functional redundancy among paralogous genes. Our previous genetic analysis in zebrafish identified important roles for HoxB- and HoxC-related genes in specifying anterior vertebral identities, whereas the contributions of HoxA- and HoxD-related genes remained unresolved. Here, we examined adult zebrafish carrying extensive combinations of hoxaa, hoxab, and hoxda cluster deletions and generated five-gene homozygous mutants carrying frameshift mutations in hoxa3a, hoxa4a, hoxa5a, hoxd3a, and hoxd4a. X-ray micro-CT analysis revealed no obvious alterations in anterior vertebral morphology in either the compound cluster mutants or the five-gene mutants. These results indicate that HoxA/D genes make only a limited detectable contribution to anterior vertebral patterning in zebrafish. Together with our previous findings, they suggest that vertebral patterning functions are distributed unevenly among zebrafish Hox clusters, with a predominant contribution from HoxB/C-related genes.

The foundation of coral reef ecosystems centered around the nutritional relationship between corals and intracellular algal symbionts. Although these symbioses are highly obligate for coral hosts, many partnerships are re-established anew with each coral generation. Furthermore, climate change destabilizes the symbiosis, and the cellular mechanisms underlying successful symbiont colonization of hosts and host development, and how they are affected by thermal stress are poorly understood. Here, we explored the effect of algal species and thermal treatments on symbiont and host cell proliferation by offering Acropora tenuis larvae one of four algal species pre-exposed to elevated or ambient temperature. In addition, we characterized the cell-surface glycome of each species-temperature combination to understand its role in symbiont recognition and proliferation. We found that thermal pre-treatment negatively affected algal photosynthetic efficiency and initial symbiont density in hosts, but did not affect symbiont colonization rate or cell proliferation. In contrast, host cell proliferation was affected in a species-specific manner. Thermal pre-treatment of B. minutum and D. trenchii negatively affected host cell proliferation compared to control symbionts, whereas thermal treatment of S. microadriaticum did not affect developmental outcomes. Further, uptake of thermally pre-treated D. trenchii decreased host cell proliferation below that of larvae not offered any symbionts, indicating that this relationship is costly to host development despite the high thermal tolerance of this species. Algal surface glycan composition varied across species but not by thermal pre-treatment, suggesting reductions in density of thermally pre-treated algae may be due to changes in physiology rather than altered surface chemistry. Further, variation in glycan abundance across species did not track differences in colonization rate or symbiont density, hinting towards a smaller role of glycans in host-symbiont specificity.

The production of gametes is one of the universal processes of life and given millennia of evolution many of the core genes involved are expected to be highly refined and resistant to change. Here we examine the pattern of expression and functionally characterize three genes widely involved in gametogenesis: Vasa, Piwi, and Dnmt1. We examine their cellular location during both oogenesis and spermatogenesis in Oncopeltus fasciatus, a hemipteran insect, using fluorescent in situ hybridization chain reaction. In females, expression of Piwi was restricted to the trophocytes, but Vasa and Dnmt1 were expressed in developing ooctyes as well as the trophoctyes. In males, Piwi was expressed in the testis germline stem cells (GSCs), but Vasa and Dnmt1 expression was absent from these cells. Additionally, Piwi and Vasa were expressed in secondary spermatogonia and spermatocytes while Dnmt1 was predominantly expressed in the primary spermatogonia. We also functionally characterized their necessity for gametogenesis after RNAi-mediated gene expression knockdown. Vasa was not required for oogenesis but was required during spermatogenesis. Piwi and Dnmt1 were required for both oogenesis and spermatogenesis. These results were in contrast to their requirement for these processes from other organisms, which highlight the unexpected variation found in the processes and suggest that conservation depends on the level at which genes are examined: sequence, expression, or phenotypic and biochemical function. This suggests there is a more nuanced evolutionary story of conserved, yet plastic, gametogenic gene set in Metazoa that deserves further study across more organisms.

Background and Aims: Dioecy has evolved independently multiple times in Palms, yet the genetic mechanisms underlying sex determination remain poorly understood across the family. Current knowledge is limited to genera from the subfamily Coryphoideae, in which Phoenix and Kerriodoxa independently recruited the same genomic region as the sex determining locus, suggesting a potentially convergent evolutionary pathway. Whether this pattern extends across the family remains unknown. Here, we investigate the genetic basis of sex determination in Phytelephas aequatorialis, from the Ceroxyloideaea palm subfamily, to identify its sex chromosome system, sex-linked genomic region, and candidate sex-determining genes. Methods: Using exon capture sequencing data from male and female individuals of P. aequatorialis, we analyzed genome wide sequence polymorphisms to detect signatures of sex linkage. Key Results: Our analyses reveal an XY sex chromosome system in Phytelephas, similar to that reported in the few dioecious palm species investigated to date. However, sex linkage maps to a distinct genomic region, indicating an independent origin of the sex-determining locus. We further identified 16 high confidence sex-linked genes with key roles in developmental processes and hormone signaling pathways, making them key candidates for involvement in sex determination. Sex-specific genetic markers were successfully developed from a subset of these genes. Conclusions: Our study suggests that in P. aequatorialis, the independent evolution of sex separation has relied on distinct genomic regions and genes, highlighting multiple evolutionary pathways leading to dioecy across the palm family. Comparative genomic analyses across additional lineages will be essential to reconstruct the evolutionary history of sex chromosomes and sex separation in this family.

Changes in gene regulation are central to developmental evolution, yet how transcription-initiation architecture contributes to developmental expression divergence remains poorly understood. Here, we integrated chromosome-level genome assemblies, developmental RNA-seq, and STRIPE-seq-based transcription start site (TSS) profiling to examine regulatory divergence between two lineages of the Daphnia pulex complex, KAP4 and CON21, across embryonic, juvenile, and adult stages. The two lineages shared a broadly conserved developmental transcriptional program, but KAP4 consistently expressed more protein-coding genes and long non-coding RNAs than CON21. Genes with stage-specific expression were weakly expressed and showed elevated nucleotide diversity and sequence divergence relative to persistently expressed genes, consistent with relaxed purifying selection. Comparative transcriptomic analyses revealed strong stage-dependent expression divergence between the two lineages: the greatest similarity occurred during late embryonic to early juvenile development rather than mid-embryogenesis, providing limited support for a canonical developmental hourglass pattern. STRIPE-seq uncovered a conserved two-peak TSS organization upstream of annotated translation-start sites, but the relative use of distal and proximal initiation regions differed between lineages. Motif analyses further indicated that these initiation regions are associated with distinct cis-regulatory sequence environments. Promoter shape was closely linked to transcriptional output, with broad promoters showing the highest expression, peaked promoters the lowest, and intermediate promoters occupying an intermediate state. Together, these findings suggest that developmental regulatory divergence in the D. pulex complex is shaped by fine-scale remodeling of TSS usage and promoter architecture.

Phenotypic plasticity exemplifies how environmental signals can shape organismal development, yet the cellular mechanisms translating ecological cues into adaptive morphologies remain incompletely characterized. Many species of Daphnia (freshwater crustaceans) develop a diverse array of inducible defenses in response to predator chemical cues (kairomones), providing a tractable system for examining how ecological pressures are transduced through developmental mechanisms. Daphnia lumholtzi produce elongated head and tailspines in response to kairomones. How this occurs at the cellular level, through changes in cell size, proliferation, or both, is currently unknown. To address this question, we quantified the temporal dynamics of cell proliferation in D. lumholtzi across 72 hours following kairomone exposure using EdU incorporation (marking proliferating cells) and DAPI staining (to quantify total nuclei). Predator cue exposure induced a three-phase proliferative response: initiation within 24 hours (slightly increased proliferation and total cells), a transitional plateau at 48 hours (minimal effects in both treatments), and commitment by 72 hours (strong increases in both proliferation and cell accumulation). Headspines exhibited higher proliferation than tailspines, suggesting anterior- posterior developmental prioritization. Treated animals maintained smaller average cell sizes throughout the response, consistent with continuous addition of newly divided cells rather than cell enlargement. Unlike the delayed-division strategy in D. longicephala or bilayer formation in D. pulex, D. lumholtzi employs sustained hyperplasia, demonstrating that the same selective pressure produces similar ecological outcomes through mechanistically distinct developmental programs. Our findings bridge ecological signals with cellular responses, exemplifying eco-evo- devo integration whereby environmental pressures (predation) are transduced through developmental processes (cell proliferation dynamics) to generate adaptive and heritable morphological diversity across ecological and evolutionary timescales.

Specialized suction discs are functional innovations that enable fishes to attach to diverse surfaces and resist hydrodynamic forces. The adhesive discs of clingfishes (Gobiesocidae) vary in size and shape, but it is unclear how ecological factors have influenced their morphological evolution. Here we analyzed the disc and skull shape of 74 clingfish species using micro-CT scanning and 3D geometric morphometrics to investigate the role of habitat and substrate use on patterns of diversification. We also present novel comparisons of adhesive performance for 10 clingfish species. Clingfish interface directly with their environment using their adhesive discs, but we found that the disc and skulls are evolutionary integrated and share similar responses to habitat. Transitions from coastal habitats to coral reefs promoted elevated rates of evolution and morphological disparity across the body, whereas transitions to freshwater did not. Concurrently, repeated shifts from living on hard substrates to softer substrates (i.e., macroalgae and seagrass) were associated with more constrained disc shapes, convergent morphologies, and differences in adhesive performance. We propose that habitat and substrate use make complementary contributions to skeletal diversification in clingfishes, but the adhesive system requires further investigation to disentangle the complex interactions between form and function.

The highly conserved body plan of nematodes makes members of this phylum excellent models to study cell type evolution. Early branching nematode lineages, mostly occupying aquatic habitats, usually contain caudal glands deployed for underwater attachment to a substrate, but have been thought to lack phasmid sensory organs, resulting in their historical classification as "Aphasmidia". With the transition to a terrestrial environment, nematodes lost caudal glands and gained phasmid sensory neurons. The supposed mutually exclusive existence of caudal glands and phasmids has led to the suggestion that phasmid neurons may have evolved from caudal glands. Here, we rule out this possibility through light and electron microscopical analysis of Mononchus aquaticus, a member of the early branching Dorylaimia lineage, showing that phasmid sensory neurons and caudal glands do coexist. This observation not only argues against a proposed cell type evolution scenario accompanying aquatic-to-terrestrial transitions but also indicates that the presence of phasmid sensory organs may have been an ancestral trait of the nematode phylum.

One of the first genes identified to regulate the spacing of nuclei in the multinucleated myofiber was Kinesin-1. However, the mechanism by which Kinesin-1 or other kinesins regulate myonuclear spacing is not known. Critically, the myofiber lacks centrosomes, and the many myonuclei act as the primary microtubule organizing centers of the cell. Because of this unique re-structuring, we hypothesized that the kinesins that drive centrosomes apart during mitotic spindle elongation may play a similar role in spacing myonuclei. We found that the bipolar Kinesin-5 (Klp61f) and the (-)-end directed Kinesin-14 (ncd) were both necessary for myonuclear spacing at different times, with both being necessary during embryogenesis, but only ncd being necessary in the fully differentiated myofiber. To investigate the shared mechanisms during embryogenesis, we used live-imaging and found that, similar to the mitotic spindle, Klp61f acts as an accelerator for myonuclear movement, whereas ncd acts as a brake contrary to this movement. To investigate these mechanisms and test the hypothesis that this is dependent on microtubule-sliding, we used super-resolution microscopy to visualize and quantify the microtubule network in animals with disrupted Klp61f or ncd. We found that in both cases, there was a decrease in the amount of microtubule overlap between neighboring myonuclei. Furthermore, we found that disrupting ncd led to severe changes in microtubule network organization, supporting our hypotheses that microtubule-sliding is necessary to space myonuclei, and that ncd likely functions through a unique mechanism in the differentiated myofiber to maintain myonuclear spacing. Together, our data supports a model where myonuclear spacing is regulated by a counteracting force generated by different kinesins during embryonic development. Furthermore, one kinesin, ncd, is repurposed in the differentiated myofiber to dynamically crosslink microtubules, a function necessary to anchor nuclei in place. Thus, kinesin motors regulate myonuclear spacing across developmental time by leveraging opposing forces through diverse mechanisms.

Stump sprouts may originate from pre-existing preventitious buds or adventitious buds, but whether these developmental origins produce distinguishable whole-shoot phenotypes remains unclear. We examined all visible sprouts on a single stump of an ornamental cherry (Prunus sp.) and measured four preselected traits: fresh weight, basal diameter, shoot length, and node number. The initial census comprised 35 sprouts. These traits were strongly correlated and were dominated by a common multivariate growth axis. Unsupervised mixture modelling identified statistical components within the morphological distribution, but these components could not be assigned to preventitious or adventitious developmental origins. During subsequent examination, we identified an unexpected basal longitudinal cavity and recorded its presence and axial length. The cavity occurred in 23 of 31 sectioned sprouts (74%). After complete removal of the initial sprout population, all visible sprouts were collected again seven days later. This second cohort comprised 22 substantially smaller sprouts but reproduced the coordinated whole-shoot growth structure. Remarkably, the basal longitudinal cavity occurred in 16 of 22 sprouts (73%), nearly identical to its initial frequency despite the marked difference in shoot size. Cavity presence was largely independent of overall shoot size, whereas cavity length increased with size among cavity-positive sprouts. Thus, gross quantitative morphology did not discriminate the hypothesized developmental origins, whereas basal anatomy revealed a reproducible feature potentially informative of developmental history. Direct anatomical and ontogenetic analyses are required to determine whether the basal longitudinal cavity is associated with sprout origin.

Tessellated calcified cartilage (TCC) is a hallmark of the chondrichthyan skeleton, yet its development early in ontogeny across the four major groups (batoids, galeomorphs, squalomorphs, and holocephalans) remains poorly understood. Specialised traits of TCC, such as multi-layered TCC and internal mineralised trabeculae, typically develop in response to feeding mechanics. In this study, we evaluated TCC morphology in the jaws of 12 representative taxa to observe its structure at an early ontogenetic stage to determine whether these specialised features had yet developed. Batoids consistently exhibited well-developed, homogeneous, polygonal tesserae early in ontogeny regardless of jaw morphology or feeding habit. In contrast, galeomorphs displayed high morphological heterogeneity. Notably, we document the first report of an extensive internal trabecular network in a non-batoid elasmobranch, observed in Ginglymostoma cirratum, which may serve to resist the mechanical pressures of specialised suction feeding. Furthermore, we identified voussoir tesserae in galeomorphs for the first time, extending their documented presence across all elasmobranch groups, where they display an inverted aspect ratio (wider than tall) compared to mature forms. The durophagous Mustelus mustelus exhibited surprisingly poor TCC development despite being a durophagous feeder, pointing to a pronounced ontogenetic lag. In Squatina oculata, TCC was characterised by large and thick tesserae and extensive fused tesseral regions which may relate to its explosive ambush predation mode, whereas the holocephalan Chimaera exhibited a poorly mineralized, mesh-like structure without resolvable discrete tesserae or trabeculae-matching findings from previous studies. Across all specimens, multi-layered TCC was absent, confirming that multi-layering develops later in ontogeny. These results demonstrate that generalised models of TCC development based on one group or a few taxa fail to capture the broader diversity of TCC morphology. It also opens up many exciting avenues for further study, and forms the basis for comparisons with fossil chondrichthyans, to investigate the evolution of TCC.

An animal's musculature must support its specific biomechanical needs, so muscle morphology and volume allocation may adapt when locomotor strategies diversify. We examined muscle size and morphology in two sister families of bombycoid moths, wild silkmoths (Saturniidae) and hawkmoths (Sphingidae), that have diverged in wingbeat frequency, wing morphology, and behavior. Although both families rely on the same muscles to power and steer flight, they may distribute muscle volume differently to prioritize distinct functions. We hypothesized that flight power muscle proportions are larger in hawkmoths and increase with wingbeat frequency, helping meet inertial power demands of high-frequency maneuverable flight. We also hypothesized that some individual muscles diverge in proportional volume and area to support distinct wing control strategies. To test our hypotheses, we took CT scans of twenty bombycoid species and quantified volumes and geometries of six flight muscle pairs. As expected, flight power muscle proportions positively correlate with wingbeat frequency and are generally greater in hawkmoths. Two of three steering muscles diverge substantially in relative volume and area between families. Most muscles exhibit greater length in silkmoths and greater cross-sectional area in hawkmoths. Finally, the dorsal oblique(DO) muscle diverges exceptionally in size and morphology, being highly developed in hawkmoths and smaller or absent in silkmoths. This unexpected difference supports the DO having an underappreciated role in flight control, possibly via shaping indirect strain propagation in the elastic thorax. We show that muscle volume distribution parallels bombycoids' divergent flight strategies, demonstrating how muscle allocation can adapt for specialized functional goals.

Nesidiocoris tenuis is an important zoophytophagous mirid bug used as a biological control agent in agriculture, and breeding efforts based on genomic information aim to increase its utility. Visible eye-color mutants are useful genetic markers because they are easily distinguishable and are therefore widely used in insect genetics and genome editing studies. Here, we investigated the genetic basis of a spontaneous red-eye mutant identified in a laboratory strain of N. tenuis. Classical crossing experiments suggested that the red-eye phenotype is controlled by a single recessive locus. RNA-seq and RNA interference (RNAi) analyses identified scarlet and cinnabar as the primary candidate genes associated with the phenotype. Further genomic analysis revealed a large deletion and insertion within exon 5 of the mutant scarlet allele, potentially causing exon skipping and disrupting transporter structure. The insertion pattern is consistent with a microhomology-mediated break-induced replication (MMBIR)/fork stalling and template switching (FoSTeS)-like event that may have been generated through polymerase{theta} -mediated repair. Together, these findings identify the causative mutation underlying the red-eye phenotype and provide a useful visible marker for future functional genetic studies and genome-assisted breeding in N. tenuis.